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We propose a mean-field model to explore the equilibrium coupling between micelle aggregation 
and lattice choice in neutral copolymer and selective solvent mixtures. We find both thermotropic 
and lyotropic transitions from face-centered cubic to body-centered cubic ordered phases of spher- 
ical micelles, in agreement with experimental observations of these systems over a broad range of 
conditions. Stability of the non-closed packed phase can be attributed to two physical mechanisms: 
the large entropy of lattice phonons near crystal melting and the preference of the inter-micelle re- 
pulsions for the body-centered cubic structure when the lattice becomes sufficiently dense at higher 
solution concentrations. Both mechanisms are controlled by the decrease of micelle aggregation 
and subsequent increase of lattice density as solvent selectivity is reduced. These results shed new 
light on the relationship between micelle structure - "crewcut" or "hairy" - and long-range order 
in micelle solutions. 

PACS numbers: 



I. INTRODUCTION 

Well known for their ability to self-assemble into a 
seemingly endless array of nanoscale structures 0, 
neutral block copolymer and solvent mixtures have re- 
cently received significant attention as a test bed for phe- 
nomena of a more classical nature: the crystallization of 
spherical particles [1 1 1, S 0, B 1 EH- In se- 
lective solvent copolymer chains aggregate into spherical 
micelles - among other structural motifs - protecting the 
insoluble polymer blocks within a core domain while al- 
lowing the soluble blocks to dissolve in a diffuse coronal 
domain. At low concentrations when inter-micelle con- 
tact is rare, correlations are fluid-like. However, at suffi- 
ciently high concentrations the swollen coronal brushes of 
neighboring micelles come into contact, and to minimize 
the repulsive effect of this contact, the micelles adopt 
lattice configurations. Based on a series of experimental 
studies of lattice-forming copolymer micelle solutions, a 
basic thermodynamic picture has emerged: generically, 
so-called "crew cut" micelles characterized by short coro- 
nal brushes favor the face-centered cubic (fee) structure, 
while "hairy" micelles with larger coronae favor the body- 
centered cubic (bec) lattice arrangement 0, 0, [1] . 

At low temperatures, many repulsive particle systems 
favor close-packed lattice structures - such as fee - pri- 
marily due to the more distant nearest neighbors and sub- 
sequent weaker repulsions in these lattices. As tempera- 
ture is raised, the entropic contribution to the free energy 
from lattice fluctuations eventually becomes comparable 
to repulsive contributions. Due to a relative abundance 
of soft-phonon modes of bee lattice, this structure is par- 
ticularly stable near to the crystal melting transition [12j • 
Indeed, quite generic consideration suggests that the bec 
structure is always preferred near to the lattice melting 
transition due to its large entropy (l3j . Entropic consid- 
erations, however, are not necessary to drive a fcc-to-bcc 
structural transition in systems of purely repulsive parti- 
cles. When repulsions are characterized by some length 



scale - as in Yukawa systems - energetic considerations 
alone show that the generic preference for close-packed 
structures at low particle densities gives way to a pref- 
erence for bec at hig h densities which persists even at 
zero temperature 0, EH, EH- Indeed, the bec lattice 
is well-known to be the three-dimensional ground state 
for Coulomb repulsions, when interactions exist over an 
infinite range [l?], EH • 

In a series of systematic studies, Lodge and coworkers 
identified fcc-to-bcc transitions in copolymer micelle so- 
lutions which are both thermotropic [a, lal _ driven by in- 
creasing temperature, T - and lyotropic\Q, @] _ driven by 
increasing polymer concentration, <i>. The former tran- 
sition scenario might suggest an entropic stabilization of 
the bec structures due to lattice fluctuations, while the 
latter scenario seems to implicate a density-driven tran- 
sition in which positional entropy does not play a crucial 
role. In this article, we construct a mean-field model to 
probe the relationship between the thermodynamics of 
micelle aggregation and the physical mechanisms under- 
lying structural phase transitions of micelles in copolymer 
solutions at fixed concentration and temperature. While 
simulations of ordered structures in star-polymer solu- 
tions [19l . |2fJ performed at fixed aggregation number, p, 
and number density, n, shed some light onto the related 
problem of copolymer micelles in solutions, the thermo- 
dynamics of the present system has the complication that 
intra-micelle energetics must be in equilibrium with inter- 
micelle effects. That is, in copolymer solutions as $ and 
T are varied, the system simultaneously adjusts the mean 
values of n and p as well as the structural organization 
of micelles. 

To address this issue, the present model combines 
treatments of good and bad solvent polymer domains 
with an effective theory for the free-energy cost - from 
micelle interactions and lattice entropy - of confining mi- 
celles within a lattice. Therefore, we address, at a mean- 
field level, how the thermodynamics of the internal de- 
grees of freedom of the constituent "particles" of this 
system, the micelles, drive the changes in inter-micelle 
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order observed on larger length scales. The principle 
consequence of increasing temperature in these systems 
is to decrease the solvent preference for coronal polymer 
domain leading to a decrease in the number of chains 
per aggregated micelle, p Q. We argue below that the 
dissolution of micelles with increased T has two effects: 
to reduce the strength of the effective micelle repulsions 
and to increase the packing density of micelles. Both ef- 
fects tend to enhance the stability of the bec structure as 
solvent selectivity is decreased, although the two mech- 
anisms are of qualitatively different origins. In the first 
case, occurring at low lattice entropy stabilizes the bec 
phase for sufficiently weak micelle interactions, while in 
the second case, at higher concentration, it is the micelle 
potential itself which prefers the bec structure over fee at 
high micelle densities. Varying the relative lengths of the 
core and coronal polymer blocks within the model, we in- 
vestigate the stability of ordered fee and bec phases as a 
function of the chain architecture. While the "crewcut" 
vs. "hairy" micelle distinction roughly describes ordered 
phase behavior, we find that the relative sizes of corona 
and core do not wholly determine lattice choice. Both 
thermotropic and lyotropic fcc-to-bcc transitions are pre- 
dicted over a broad range of polymer architectures. 

This article is organized as follows. In the Sec. UT1 we 
construct a mean-field model for the free-energy of the 
internal state of the micelles, as well as the cost of lattice 
confinement. We present results for this model in Sec. 
IIIII and conclude with a brief discussion in Sec. IIVI 



II. FREE ENERGY MODEL FOR LATTICE 
MICELLE PHASES 

In this section we construct an effective mean-field 
model to compute the free energy of lattice micelle phases 
of copolymer solutions at fixed temperature and polymer 
fraction, $. We first consider the energetics of isolated 
micelles in solution before turning to the inter-micelle ef- 
fects. Note that the model for isolated micelle energetics 
presented here has been recently analyzed in detail in the 
limit of low concentration (in the absence of inter-micelle 
effects) by Zhulina et al. [H]]. We present a complete 
discussion here for the purpose of elucidating the results 
of the full calculation, which takes into account effects of 
lattice confinement. 



A. Internal Micelle Configurations 

We consider a solution of diblock copolymers of N to- 
tal monomers of volume fraction <&. The lengths of the 
insoluble core and soluble coronal blocks are 7V C and Nt, 
respectively. Temperature dependence enters the model 
through the selectivity of the solvent for the coronal 
brush domain over the core domain. Here, we assume 
that solvent is excluded from the core domain entirely, so 
that the insoluble polymer blocks are exposed to unfavor- 



able solvent contact only at the interface separating the 
core and coronal domains. Microscopic solvent-monomer 
interactions give rise to a surface tension, 7, for the core- 
corona interface. We generically expect the importance 
of molecular scale interactions to decrease with temper- 
ature, and hence, 7 to be a decreasing function of T as 
is typically the case for example, in immiscible polymer 
blends [22J . In this model the effect of solvent selectivity 
will be encoded entirely in 7. The free energy of an iso- 
lated micelle configuration, Fq, is the sum of three parts, 

F = F c + F mt + F b , (1) 

where F c and Ft, are the free energies of the core and 
coronal brush domains and Fi nt is the free energy cost 
of exposing the core domain to solvent along the inter- 
domain surface. 

Consider an isolated micelle composed of p total chains 
shown in Figure [TJ The core domain extends to a radius, 
R c , while the tips of the coronal brush reach to a radial 
distance, i?o- Because solvent is absent from the core, the 
polymer blocks that constitute this domain are described 
by the statistics of the molten polymer brush [23| . As 
such, the free energy of this domain can be calculated in 
the strong-stretching approximation for spherical, molten 
polymer brushes [24j . 

_F^ = 3fp_Rl_ 

k B T 80 N c a 2 ' 1 ' 

where a is the mean monomer-monomer distance. 
Since excluded-volume interactions between chains are 
screened in the molten polymer state, the free energy 
of this domain is purely due to the entropic (Gaussian) 
cost of extending the chains in the brush. Given the sur- 
face tension between core and coronal brush domains, we 




FIG. 1: A cartoon of an isolated micelle in solution. Solvent 
is excluded from a molten polymer core of radius, R c . A 
dissolve coronal brush extends from R c to a radius of Ro- In 
the Daoud and Cotton model (see text) this brush can be 
divided into correlation "blobs" of size, £, determined by the 
distance between neighboring chains. 
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have, 



Ml 



(3) 



Treating the core block as hydrophobic, we expect that 7 
is on the order of 10 mN/m (Tlj . which in thermal units 
corresponds to 7a 2 /ksT « 0.5, assuming a segment size 
of a = 5A. Note that from the volume constraint on the 
core we have, 
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where po is the monomer density. 

The coronal domains of the micelles are under good- 
solvent conditions, leading to a dramatic swelling of these 
blocks in comparison to the core domain. We treat the 
free energy of this domain within the Daoud and Cotton 
model of "star-like," spherical polymer brushes [2511. This 
model has been studied extensively elsewhere [2a l27j ; 
we present a detailed discussion here since the physics 
of polymer corona has important implications for the 
micelle-micelle interactions. Under good-solvent condi- 
tions, the properties of a semi-dilute polymer solution are 
governed largely by a single length scale, £, the so-called 
blob size, which is determined by the mean distance be- 
tween neighboring chains (23[. At shorter length scales 
than £ the correlations of the chains obey the Flory scal- 
ing, g 3 w po^/a 2 , where g is the number of monomers 
per blob. For a spherical brush of p chains, the aver- 
age lateral distance between neighboring chains is deter- 
mined by the geometrical condition that Anr 2 = p£, 2 (r) 
(see Figure HJ so that, 



(5) 



The local concentration of monomer in the dissolved 
brush is determined from 4>b(r) — g(r)/£, 3 (r) oc p 2 ^ 3 /r 4 / 3 . 
Because the total number of monomers in the corona is 
fixed, integrating this concentration profile over the vol- 
ume of the coronal brush allows us to solve for the size 
of the brush, 
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pNb — 4-7T / dr r 2 (j)b(r) 



(6) 



Inverting this condition we find that the size of the mi- 
celle obeys, 
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(8) 



and i?„ is the size of a star polymer of p chains of length 
Ni, in the limit that R c — > 0. Note that in the "hairy" 



micelle case, when i?» 3> R c , we have Rq 
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FIG. 2: The solid line shows the crossover in the preference 
of U(r) for the bec and fee lattice structures as function of 
micelle aggregation, p, and reduced density, 77 = (n/6)na 3 . 



whereas for "crew-cut" micelles, when R* < i? c , we have 
Ro ~ (Nop) 1 / 3 . 

The free energy of the good solvent domain is com- 
puted by counting the total number of these thermal 
blobs, each of which represent a fcgT worth of energy. 
This is most conveniently accomplished by integrating 
the blob density, £,~~ 3 (r), over the coronal domain vol- 
ume, from which it can be shown, 
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In the absence of any inter-micelle effects, the scaling be- 
havior of the mean aggregation number can be obtained 
by considering the "hairy" and "crew-cut" limits of Fq 
|27| . When R a ^> R c , the free energy of the coronal 
brush, Fb, dominates over the core stretching energy, F c , 
so that p r~j 7 6 / 5 iVc^ 5 . In the opposite crew-cut limit, 
Fb <C F c , the mean aggregation grows as p ~ 7-/V c . De- 
creasing solvent selectivity (decreasing 7), decreases the 
mean number of chains per micelle, yielding a roughly lin- 
ear dependence of chain aggregation on the core-corona 
surface tension. 



B. Micelle Repulsion and Lattice Entropy 

At low polymer concentrations, the density of micelles 
is sufficiently small that inter-micelle effects may be ig- 
nored [28]. As $ is increased well beyond the critical 
micelle concentration the coronal domains of neighbor- 
ing micelles begin to overlap, leading to strong repulsive 
effects. To minimize the free-energy cost of such repul- 
sions, micelle solutions adopt periodic arrangements, al- 
ternately bee or close-packed structures. We consider a 
lattice of micelles of aggregation number, p, at number 
density, n. Given the low solubility of single, unaggre- 
gated chains in solution [H, [2i| , we assume a negligible 
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(zero) concentration of free chains. Therefore, a mean 
micelle aggregation number of p implies a micelle den- 
sity, 



pN 



(10) 



Note that as micelles dissolve and p decreases, the micelle 
density increases to maintain a fixed polymer fraction, 

When neighboring micelles begin to overlap at high 
densities (i.e. n > (2i? ) -3 )> we must account for the 
repulsive effects between coronal brushes. By a scaling 
analysis, Witten and Pincus showed that the pair poten- 
tial for star-like, spherical brushes exhibits a logarithmic 
interaction at short distances [2(|. Physically this is jus- 
tified by the demand that the free-energy cost of bring- 
ing two micelles of p chains each to a separation of order 
2R C , should be roughly the same as the free energy of 
creating a single micelle of 2p chains (see eq. ©). On 
these grounds, Likos and coworkers proposed the follow- 
ing form of the pair potential for star polymers of p chains 



the bcc or fee structures as a function of p and the re- 
duced density, r\ = (ir/6)na 3 . At low densities, n < 0.7, 
micelles in either lattice interact via the long-distance 
portion of the potential. In this regime, lattice geom- 
etry dictates that the potential energy of fee structure 
is exponentially smaller than that of the bcc structure 
due to its more distant nearest-neighbors. From studies 
of Yukawa systems [1J, [l5j, [l6| , it is known that the bcc 
structure is generically favored at high densities, and this 
fact accounts for the p-dependent crossover in preferred 
structure shown in Figure [5] for low aggregation number 
(p < 100). Near to r\ « 0.7 nearest-neighboring micelles 
reach the short-distance, logarithmic regions of U (r), and 
the crossover density becomes only weakly dependent on 
p. The crossover from preferred fee to bcc structures with 
increased density is a rather generic property of repulsive 
systems, qualitatively insensitive to the precise form of 
the inter-micelle potential. 

We compute the lattice entropy of competing struc- 
tures by considering the harmonic free-energy cost of dis- 
placing micelles by Ar" from their equilibrium positions, 
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At shorter separations than a soft-core size, ex, the repul- 
sion has the Witten and Pincus form. Since the micelles 
have a finite size, at separations beyond a, U(r) has a 
Yukawa form, which falls off exponentially on a length 
scale set by the largest blob size, £ ma x (see Fig. [J). Note 
that both the potential and force are continuous at r = a. 
The crossover in the potential is determined by the radial 
separation at which the outermost blobs of neighboring 
micelles fully overlap. From this we have, 
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This pair potential has been shown to agree well with 
scattering experiments in dilute solutions of copolymer 
micelles [10|, EJ and star polymers [3(| • 

To compute the free-energy cost of coronal overlap 
for micelles in the various structures, we perform lat- 
tice sums of the potential in eq. (|11|) numerically. Be- 
cause U (r) falls off exponentially with distance, the sums 
are truncated to include a finite number of lattice neigh- 
bors, yielding a numerical precision of better than one 
part in 10 s . Unlike the case of pure Yukawa systems 
which are characterized by a single length scale, the 
inter-micelle potential depends explicitly on two separate 
length scales, namely, a and £ max . To analyze the influ- 
ence of the inter-micelle potential on equilibrium struc- 
ture, in Figure Owe plot the preference of U(r) for either 
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where the i and j label Cartesian directions, a and (3 label 
lattice sites and DfJ 3 is the dynamical matrix computed 
from, 
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with d/drf the derivative with respect to the i compo- 
nent of the position of the ath lattice site. This matrix 
can be be diagonalized in Fourier space to find the eigen- 
values, Aj(k), associated with the three phonon polar- 
izations at wavevector k. In the harmonic approxima- 
tion [3~H the equilibrium mean-square displacement can 
be computed by 
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where {G} is the set of reciprocal lattice vectors corre- 
sponding to a lattice of N m — nV micelles. The free- 
energy contribution (per micelle) from positional fluctu- 
ations in the harmonic approximation is computed by, 
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Although the potential in eq. (fTl"j) and its first deriva- 
tive are continuous over the full range, the second deriva- 
tive is discontinuous at r = a, which would enter the 
computed phonon spectrum through eq. (|15p . Hence, 
the form of U{r) proposed in ref. [30j, would lead to 
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certain anomalies in the dependence of the fluctuation 
contribution to the free energy on micelle density, n. To 
circumvent this issue while maintaining some of the de- 
tailed dependence of the vibrational spectrum on the 
inter-micelle potential, we use only the long-distance, 
Yukawa form of U(r) - even down to r < a - when 
computing the dynamical matrix. We expect that this 
becomes a poor approximation well above the density at 
which the nearest-neighbor separation reaches a. Hence, 
in this model the vibrational properties of the micelle lat- 
tices will be identical to those characterized by Robbins, 
Kremer and Grest for lattice phases of repulsive Yukawa 
systems [TBI [l6|. 

To calculate the dependence (|Ar| 2 ) and Ffi uct on p 
and n, we compute Dfj 3 using eq. ([15)) and the long- 
distance form of U(r), numerically performing lattice 
sums. The eigenvalues A^(k) are computed, and then 
the Fourier sums in eqs. (flB]) and ([!?)) are carried out 
by summing over a finite set of reciprocal lattice vectors. 
More than 300,000 reciprocal lattice points are used to 
achieve a numerical convergence of (|Ar| 2 ) and Ffi uct to 
better than one part in 10 4 . 

When micelle interactions are sufficiently weak, lat- 
tice fluctuations will lead to a melting of the ordered 
phases. To assess this instability to the fluid micelle 
phase, we apply the Lindemann criterion for crystal melt- 
ing, which defines the melting of lattice order to occur 
when dAr] 2 ) 1 / 2 /r nn > a m , where R nn is the nearest- 
neighbor lattice spacing and a m is a numerical constant 
less than 1 |32j |. Here, we use a Lindemann parame- 
ter of a = 0.19, which was found to yield good agree- 
ment with the Hansen- Verlet criterion for meltin g |33f| in 
numerical studies of repulsive Yukawa systems [16J. As 
was found in these studies, the soft-phonon modes of the 
bcc crystal in the harmonic model are completely unsta- 
ble to thermal fluctuations for sufficiently low densities 
(i.e. n^ax ~ 0.002), and therefore, the model discussed 
above cannot be used to compute the vibrational prop- 
erties of the bcc structures at these densities. In this 
regime, we use an extrapolation of Ffi uct from the sta- 
ble bcc phase at higher densities to model the entropic 
contributions to the free energy of the lower density fluid 
phase. 

III. MODEL RESULTS AND MICELLE 
SOLUTION THERMODYNAMICS 

In this section we discuss the mean-field results of the 
model for copolymer micelle solutions introduced in Sec. 
HT1 Calculations are performed for fixed polymer content, 
and polymer composition, N c and JV&. Temperature 
dependence of the solvent selectivity is encoded in the 
core-corona surface tension, 7. For simplicity we assume 

— 1/3 

that a = p . The total free energy (per micelle), F tot 
is the sum of internal micelle free-energy contributions 
- eqs. ©, (O and §§§ - along with the inter-micelle 
free-energy contributions, U and Ffi uct . For fixed so- 
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FIG. 3: Equilibrium phase behavior for solution of symmetric 
copolymers N c = Nb = 50. The regions of stable fee and bcc 
phases are indicated, as well as the region where the crystal 
structures are unstable to the micelle fluid phase as deter- 
mined by the Lindemann criterion (see text). All phase tran- 
sitions are first order. Portions of the fcc-bcc phase boundary 
are labeled according to the physical mechanism which dom- 
inates the thermodynamics as discussed in the text. 

lution conditions, the intensive free energy per chain, 
Ftot/ Pi is minimized over micelle aggregation number, 
p. Mean-field free energy calculations for ordered fee and 
bcc micelle structures are compared to construct phase 
diagrams. We begin with a detailed discussion of micelle 
thermodynamics for solutions of symmetric copolymer 
chains. 



A. Symmetric Chains 

In Figure [3] we show the predicted phase behavior 
of copolymer solutions for symmetric molecules with 
Y c = N\j = 50. At sufficiently large core-corona surface 
tension, 7a 2 /ksT > 0.08, and polymer volume fraction, 
<f> > 0.05, the micelle solution orders into either the fee or 
bcc structures. The former structure is preferred at lower 
$, while the stable bcc structure is predicted at higher 
polymer fractions and at lower 7, before crystal melting. 
The topology of the mean-field phase diagram for the 
ordered micelle phases is consistent with experimental 
observations of symmetric diblock solutions, excepting 
here that the temperature dependence of solvent quality 
is encoded in 7 so that direct comparison to experiment 
is not possible at present. Nonetheless, the predicted 
fee— >bcc^fluid micelle phase sequence with increasing 
temperature is a generic feature observed for a variety 
of neutral copolymer solutions 0, 0, 0]- Moreover, we 
note that many systems (though not all) exhibit the fcc- 
bcc structural transition with increasing concentration as 
predicted by the mean- field model (see. e.g. ref. Q). 

To put the thermodynamics of lattice choice in the con- 
text of equilibrium micelle structure we plot the depen- 
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FIG. 4: Mean-field values of aggregation number, p, (solid 
line) and reduced density, n = (-k/6)iio~ 3 for solutions of sym- 
metric chains (N c — Nb = 50) of polymer fraction $ = 0.175. 
The arrow at ■ya 2 /ksT — 0.222 indicates the location of the 
bcc-fcc phase boundary. Hence, for values of surface tension 
smaller (larger) than this the plotted quantities correspond to 
the equilibrium bcc (fee) structure. 



dence of mean aggregation and reduced density as func- 
tions of the core-corona surface tension in FigureUJ Since 
7 is assumed to be a decreasing function of T, this phase- 
space trajectory describes a thermotropic fcc-bcc transi- 
tion. Here, we see that p has a roughly linear dependence 
on 7, not unlike the case where inter-micelle effects are 



ignored [27j. In the range 7a 2 /fcsT > 0.12, 77 is a mono- 



tonically decreasing function of 7. This stems from the 
fact that micelle density, n, is inversely proportional to p, 
while at large aggregation numbers <r 3 ~ Rq grows slower 
than linearly with increased p (for example, i?g ~ p 3 / 5 
for hairy micelles). Hence, we have the generic property 
that at large aggregation numbers increasing core-corona 
repulsion makes the system less dense on the scale of the 
intermicelle repulsion, a. Indeed, the rise in density with 
increased temperature is observed in scattering experi- 
ments from micelle lattices in solution Q . From eq. (Tl2|) 
we have a ~ RoP 1 ^ 2 for small p, thus r] decreases to 
as p — > 0, which accounts for the decrease of the reduced 
micelle density as 7a 2 /fc^T is decreased below 0.12. For 
a broad range of conditions in this model, when p is re- 
duced below roughly 20 chains per micelle, ordered struc- 
tures are unstable to the fluid phase, in agreement with 
the order of magnitude predictions for the critical chain 
number for ordering in solutions of star polymers [34j |. 

In principle, both the results that the system effec- 
tively becomes more dense and the number of chains per 
micelle decreases with decreasing core-corona surface ten- 
sion point to the bcc structure as the stable micelle lattice 
at high temperatures. The former effect drives an inter- 
micelle potential preference for bcc (see Fig. [2J while 
the latter effect suggests that lattice fluctuations stabi- 
lize this phase as micelles dissolve, since the strength of 
effective micelle repulsions is proportional to p 3 / 2 . In or- 
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FIG. 5: Differences between bcc and fee the total, AF to t, 
interaction AU, and fluctuation, AFfi uc t free energies of so- 
lutions of symmetric chains (N c = Nb = 50) as a functions 
of core-corona surface tension. The solid line corresponds to 
results at $ = 0.175 while the dashed line corresponds to 
$ = 0.125. 



der to discriminate between these two distinct physical 
mechanisms for the thermotropic fcc-to-bcc transition - 
driven by either lattice potential or lattice entropy - we 
plot the difference between F to t, U and Ffi uct (per poly- 
mer chain) for the bcc and fee structures as a function 
of 7 in Figure [5l As is the case for pure Yukawa systems 
[l6j |. the relatively large entropy of the fluctuations in 
the bcc structure always favors that structure over fee, 
so that AFfi uc t < for all conditions. However, this en- 
tropic preference amounts to roughly 10 _3 fcsT per chain, 
which should be compared to differences in lattice po- 
tential, which tend to vary on the scale of 10~ 2 k B T per 
chain. 



7 




0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 
<I> 



FIG. 6: Dependence mean aggregation on $ for symmetric 
copolymer solutions at fixed 7a 2 /knT and fixed bcc structure. 



Consider first the lower polymer fraction case, where 
$ = 0.125. Over the range of 7 plotted in Figure El 
AU is always positive indicating that interactions fa- 
vor the fee structure for this polymer fraction. Below 
7a 2 //csT = 0.105, however, the bcc structure is ther- 
modynamically stable over the close-packed lattice until 
melting for ^a 2 /k B T < 0.067. In this regime of weak 
core-corona repulsion, mean-aggregation of micelles is 
sufficiently small (p < 20) that differences in Ffi uct which 
favor the bcc lattice become comparable to the prefer- 
ence of micelle interactions favoring the fee structure. 
Hence, we view this thermotropic fcc-to-bcc transition 
to be of entropic origin, resulting from the abundance of 
soft-phonon modes in the bcc crystal. 

We now turn to the thermotropic fcc-to-bcc transition 
at $ = 0.175 shown also in Figure El In contrast to 
the lower concentration case we find a high temperature 
region, 0.057 < 7a 2 /ksT < 0.179, wherein micelle inter- 
actions prefer the bcc structure over fee. Indeed, from 
Figure 0J we find that the system becomes more dense 
as 7 is lowered, so that 77 > 0.7 in this region. While 
the generic preference of Ffi uct for the bcc structure may 
further stabilize this lattice over the fee structure, the 
plot of AU for this concentration demonstrates that the 
differences in micelle interactions dominate the entropic 
contribution to the free energy over much of the range of 
the stable bcc lattice (note that the bcc lattice melts be- 
low 7a 2 /fceT = 0.077 for this polymer fraction). There- 
fore, in this higher concentration regime the structural 
transition is more appropriately attributed to micelle in- 
teractions which favor the bcc structure as the system 
becomes effectively more dense. 

Micelle aggregation is a decreasing function of tem- 
perature, which in turn weakens micelle interactions and 
increases the number density of micelles. In Figure El 
crossover from entropy-driven to interaction-driven fcc- 
to-bcc structural transitions is marked by a cusp in phase 
boundary at $ = 0.160 and -/a 2 /k B T = 0.120. For mi- 
celles with p > 100, there is effectively a critical value of 



reduced density, rj* ~ 0.7, above which the long range of 
inter-micelle repulsion leads to a stabilization of the bcc 
structure. Since n is linearly proportional to the poly- 
mer fraction, increasing $ allows the system to pass into 
the bec-preferred region, 77 > 0.7 at larger aggregation 
number, or alternatively, larger 7. This accounts for the 
^-dependence of the fcc-bcc boundary. For <I> < 0.160, 
the weakening of the micelle potential with the loss of 
chain aggregation supersedes the density-driven transi- 
tion, leading to the entropic stabilization of the bcc phase 
near to the melting point. Therefore, we find that the 
interaction-driven fcc-to-bcc transition tends to occur at 
higher $, with the entropically-driven structural transi- 
tion occurring at lower polymer concentrations. 

Having analyzed the thermotropic transition in detail, 
we now briefly consider the case of fixed 7. Clearly, the 
discussion above implies the low-$ stability of the fee 
structure and the subsequent stability of bcc at higher 
concentrations due to the micelle interactions. However, 
the predicted fluid^bcc^fcc^bcc phase sequence with 
increasing $ predicted for the range 0.105 < 7a 2 /fc^T < 
0.265, implies an unusual thermodynamic dependence of 
micelle aggregation on <I> (shown in Figure El) ■ While at 
low concentrations p will be insensitive to <E>, when mi- 
celles begin to interact strongly the system adjusts the 
micelle structure in order to reduce the free energy cost 
of micelle overlap. In the concentration regime where mi- 
celles interact via the Yukawa portion of U(r), the mean 
aggregation increases roughly linearly with <I> in order 
to maintain a density where nearest-neighbor micelles do 
not overlap strongly. Once the system reaches a den- 
sity where nearest neighbors begin to interact via the 
short-distance, logarithmic portion of U{r), this trend 
reverses and aggregation then decreases with increased 
$. This effect not only increases the density of the lat- 
tice but also reduces the strength of repulsions. The non- 
monotonic dependence of p on <3> shown in Fig. El explains 
the re-entrant bcc phase at high concentration, since this 
structure is generically preferred at high density and low 
p. This unusual result that p first increases with $ and 
then decreases, is a fundamental consequence of the fact 
that U(r) is finite-ranged and "ultra-soft" meaning that 
it should remain finite down to very small separations. 
To prevent significant overlap, the system first slightly 
increases p, reducing density. Once strong overlap be- 
comes unavoidable, it becomes preferable to decrease p, 
since the free-energy cost of fusing two micelles cannot 
be significantly larger than creating a single micelle of 2p 
chains. 



B. Copolymer Composition: "Hairy" vs. 
"Crew-Cut" Micelles 

We now consider the influence of copolymer composi- 
tion on the predictions of the mean-field model. Defining 
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FIG. 7: Mean-field phase diagrams for asymmetric N = 100 copolymer micelle solutions. The phase are labeled as in Figure 
[3] Qualitatively, the / = 0.3 corresponds to a typical phase diagram for "hairy" micelles, while the / = 0.7 represents the 
"crew-cut" micelle case. 



/ as the fraction of chain which is insoluble, 

/=^^, (18) 

we show the mean-field phase behavior for micelle solu- 
tions of corona-rich, / = 0.3, and core-rich, / = 0.7, 
copolymers in Figure [JJ Qualitatively, these two cases 
correspond to the "hairy" and "crew-cut" descriptions of 
micelle structure, respectively. More quantitatively, we 
can characterize the micelle structure by the ratio of the 
height, H b = Rq—R c , of the coronal block to the core size. 
At $ = 0.2 and ja 2 /k B T = 0.2, we have H b /R c = 2.734 
for / = 0.3 and H b /R c = 0.768 for / = 0.7. Another 
significant difference is the mean aggregation since p is 
an increasing function of N c : p = 28.5 and p — 188, 
for the "hairy" and "crew-cut" cases, respectively. These 
theoretical phases diagrams recapitulate the observations 
of Lodge and coworkers that structural order of "hairy" 
( "crewcut" ) micelles is principally of bcc (fee) type [|| . 

It is the latter distinction in micelle aggregation for 
these two cases which accounts for the broad differences 
in the mean-field phase diagrams in Figure [JJ Scaling 
analysis of Fq suggests a roughly linear dependence of 
mean-aggregration number on the length N c , which jus- 
tifies the large differences in p for the two cases. Since 
a smaller micelle aggregation implies a larger micelle 
density, for given solution conditions "hairy" micelles 
are characterized by highly overlapping coronal brushes 
while the "crew-cut" micelles tend to form less dense lat- 
tices. In the low temperature regimes for which (p > 50 
this implies that micelle interactions tend to stabilize 
the bcc lattice in the former case and fee in the latter. 
Clearly, the effect of micelle structure does not drastically 
alter the entropic stabilization of the bcc phase near to 
lattice melting. 

In Figure[5]we plot the fixed 7a 2 /fcsT = 0.2 phase be- 
havior of N = 100 copolymer solutions as a function of $ 
and /. Here again, the core-rich copolymer micelles favor 



the close-packed structure, while the corona-rich micelles 
favor bcc. This phase diagram reveals that the re-entrant 
bcc — > fee — > bcc lyotropic phase sequence is a feature 
of copolymer solutions over a broad range of chemical 
compositions. Although the "hairy" vs. "crew-cut" dis- 
tinction broadly describes ordered-phase thermodynam- 
ics of micelle solutions, we find that equilibrium structure 
cannot be strictly determined by / (or even H b /R c ). In 
fact, basic considerations of micelle interactions dictate 
that ordered-phase thermodynamics of micelles must be 
insensitive to the size of the core domain, R c . Micelle in- 
teractions govern both repulsive and entropic free-energy 
costs of lattice formation, and from eq. (fTTj) we see that 
U(r) depends only the number of chains per micelle and 
the effective micelle size, a. Certainly, this must be the 
case, given that the self-similar statistics governing the 
monomer density in the coronal brush domain are inde- 
pendent of the length of coronal chains, Nb [IH . There- 
fore, since neighboring micelles interact through overlap 
of coronal domains where the monomer distribution is 
governed purely by p and r, such interactions cannot 
encode information about the relative sizes of the core 
domain. Or put another way, at the tips of their coro- 
nal brushes micelle interactions cannot "know" how far 
away the core domain is from the outermost edge of the 
micelle. According to a mean-field description of micelle 
solutions, differences in ordered-phase thermodynamics 
between the "hairy" and "crew-cut" cases are more prop- 
erly attributed to differences in solvation properties of the 
copolymers of differing compositions, rather than the rel- 
ative size of core and coronal domains. 



IV. DISCUSSION 

In the previous sections we have presented and ana- 
lyzed a mean-field model for the thermodynamics of con- 
centrated solutions of neutral block copolymers. This 
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analysis unifies a broad range of observed phenomena 
regarding structural transitions in lattice phases of mi- 
celles. In particular, the mean-field analysis elucidates 
the physical mechanisms responsible for phase transitions 
between ordered states of micelles. In contrast to other 
classical particle systems, such as charged-stabilized col- 
loids [HI or even hard spheres [HI, modeling the ther- 
modynamics of micellar systems has the additional com- 
plication that the internal structure of the micelles them- 
selves establishes equilibrium with the inter-micellar de- 
grees of freedom. 

This fact underlies the basic phenomenology of concen- 
trated solutions of copolymer micelles. As temperature 
is raised and solvent selectivity decreases, the mean value 
of p decreases. Since the strength of inter-micelle repul- 
sions grows as p 3 / 2 , the decrease in micelle aggregation 
leads to a rise in the relative importance of the entropic 
contributions to lattice free energy resulting from phonon 
fluctuations. As t he g eneric arguments of Alexander and 
McTague dictate 13j, this effect promotes the stability 
of the bec structure near to the micelle lattice melting. 
However, we find that even far from the regimes where 
entropic considerations are important - when solvent- 
selectivity is strong and p is large - micelle solutions 
exhibit broad ranges of bec stability. In these regions 
of phase space, stability of the bec structure over the 
close-packed fee structure can be attributed to micellar 
interactions which are of sufficiently long length scale. 

While the present model effectively describes the ther- 
modynamics associated with the fcc-bcc structural tran- 
sitions in micelle solutions, the approach is limited in 
regard to other equilibrium structures. In particular, the 
star-polymer potential of eq. (fTTj) has been shown to 
stabilize non-cubic and even diamond lattice structures 
[l9l |2C| . Furthermore, studies of similar soft-molecular 
systems suggest the stability of the relatively open A15 
lattice of micelles (3(| H3] ■ Due to the abundance of ther- 
modynamic studies of copolymer solutions exhibiting the 
structural fcc-bcc transition, we focus on this particu- 
lar phenomenon, and consequently, these lattice struc- 
tures. Generally, we might expect the stability of other 
micelle lattice phases at higher polymer concentrations. 
More likely however, is the possibility that other self- 
assembled morphologies intervene at higher concentra- 
tion, such as cylindrical, lamellar or network morpholo- 
gies. Since the inter-micellar terms in the free energy are 
necessarily of the same order as the intra-micellar ones, 
at high concentrations, the system achieves a more effi- 
cient packing by a rearrangement of aggregate topology. 
Such phase behavior is well-documented in experimental 
observations of more concentrated copolymer solutions 

Other limitations of the model include the fluctuation 
effects of the micellar structure. One of these is the 
swelling of the core domain by solvent at high temper- 
ature. Neutron scattering in dilute solutions of micelles 
reveals that penetration of solvent into the core domain 
may reach as high as 30% by volume near to micelle 




FIG. 8: Mean-field phase behavior of N = 100 copolymer so- 
lutions as a function of chemical composition, /, and polymer 
content $, for fixed solvent-core repulsion, 7a 2 /ksT = 0.2. 



lattice melting [8j. While we have not included this ef- 
fect in our mean- field treatment, we do not expect it to 
change the qualitative view of ordered phase thermody- 
namics presented here. Namely, our principle conclusion 
that structural fcc-to-bcc transitions are driven by in- 
creases in micelle density will remain unchanged. Indeed, 
the swelling of the core domain will tend to counteract 
the decrease in micelle size by loss of polymer chains as 
temperature is raised. Thus, the basic prediction that 
77 largely increases as temperature rises will hold in the 
case that solvent is permitted to swell the micelle core. 

A second type of fluctuation effect absent from the 
present model are equilibrium fluctuations in micelle ag- 
gregation number, p, leading to a certain amount of poly- 
dispersity in micelle size. As in micellar phases of copoly- 
mer/homopolymer blends [38[, fluctuations in micelle ag- 
gregation above the critical micelle concentration can be 
estimated to be on the order of 10% or less. While the de- 
gree of polydispersity in the size of micelles will be even 
smaller than this, it is not clear how even this modest 
level of polydispersity should affect the relative stability 
of various ordered phases. Certainly, a broad distribu- 
tion in micelle size and aggregation number will frus- 
trate periodic order, but given that copolymer solutions 
are known to form well-ordered structures over a broad 
range of equilibrium conditions we expect this effect may 
not be so significant for this system. 

Finally, we conclude by contrasting the thermodynam- 
ics of micellar phases to other spherically-aggregating 
molecular systems, such as ionic or non-ionic surfac- 
tants. Indeed, surfactant systems exhibit structural tran- 
sitions between competing ordered structures as concen- 
tration and solvent quality are varied [39| . However, 
unlike their polymeric counterparts, aggregates of low- 
molecular weight surfactants tend be characterized by 
a fixed size, dictated by molecular-scale packing con- 
straints [40j. In contrast, micelles of flexible polymeric 
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amphiphilcs adjust their aggregation by a factor 5-10 
times over the accessible range of solution conditions. 
This molecular flexibility is fundamentally linked to the 
unusual feature that copolymer micelle lattices become 
more dense as temperature is raised and molecular forces 
are weakened, ultimately driving transitions in the long- 
ranged structural order of the system. 



supported by NSF Grants DMR01-02459 and DMR04- 
04507, the Donors of the Petroleum Research Fund, Ad- 
ministered by the American Chemical Society and a gift 
from L.J. Bernstein. 



Acknowledgments 

It is a pleasure to thank J. Bang, R. Kamicn, P. Zihcrl 
and C. Santangelo for helpful discussions. This work was 



[20] 



D. E. Discher and A. Eisenberg, Science 297, 967 (2002). [21 
S. Jain and F. S. Bates, Science 300, 460 (2003). 
G. A. McConnell, A. P. Cast, J. S. Huang and S. D. [22 
Smith, Phys. Rev. Lett. 71, 2102 (1993). 

G. A. McConnell and A. P. Cast, Phys. Rev. E 54, 5447 [23 
(1996). 

J. Bang,T. P. Lodge, X. Wang, K. L. Brinker and W. R. [24 
Burghardt, Phys. Rev. Lett. 89, 215505 (2002). [25 
T. P. Lodge, B. Pudil and K. J. Hanley, Macromolecules 
35, 4707 (2002). [26 
J. Bang, K. Viswanathan, T. P. Lodge, M. J. Park and 
K. Char, J. Chem. Phys. 121, 11489 (2004). [27 
T. P. Lodge, J. Bang, M. J. Park and K. Char, Phys. 
Rev. Lett. 92, 145501 (2004). [28 
T. P. Lodge and J. Bang, Phys. Rev. Lett. 93, 245701 
(2004). [29 
J. Stellbrink, G. Rother, M. Laurati, R. Lund, L. Willner [30 
and D. Richter, J. Phys.: Condens. Matter 16, S3821 
(2004). 

M. Laurati, J. Stellbrink, R. Lund, L. Willner, D. Richter [31 

and E. Zaccarelti, Phys. Rev. Lett. 94, 195504 (2005). 

W. G. Hoover, D. A. Young, and R. Grover, J. Chem. [32 

Phys. 56, 2207 (1972). [33 

S. Alexander and J. McTague, Phys. Rev. Lett., 41, 702 [34 

(1978). 

D. Hone, S. Alexander, P. M. Chaikin and P. Pincus, J. [35 
Chem. Phys. 79, 1474 (1983). 

K. Kremer, M. O. Robbins and G. S. Grest, Phys. Rev. [36 
Lett. 57, 2694 (1986). 

M. O. Robbins, K. Kremer and G. S. Grest, J. Chem. [37 
Phys. 88, 3286 (1988). [38 
K. Fuchs, Proc. Roy. Soc. (London) A151, 585 (1935). 
W. J. Carr, Phys. Rev. 122, 1437 (1965). [39 
M. Watzlawek, C. N. Likos and H. Lowen, Phys. Rev. 
Lett. 80, 4450 (1999). [40 
C. N. Likos, N. Hoffmann, H. Lowen and A. A. Louis, J. 
Phys.: Condens. Matter 14, 7681 (2002). 



E. B. Zhulina, M. Adams, I. LaRue, S. S. Sheiko and M. 
Rubinstein, Macromolecules 38, 5351 (2005). 

F. S. Bates and G. H. Fredrickson, Annu. Rev. Phys. 
Chem. 41, 525 (1990). 

P.-G. de Gennes, Scaling Concepts in Polymer Physics, 

(Cornell University Press: Ithaca, 1979). 

A. N. Semenov, Sov. Phys.-JETP 61, 733 (1985). 

M. Daoud and J. P. Cotton, J. Phys. (Paris) 43, 531 

(1982). 

T. A. Witten and P. A. Pincus, Macromolecules 19, 2509 
(1986). 

T. M. Birshtein and E. B. Zhulina, Polymer 30, 170 
(1989). 

A. Johner, J. F. Joanny and C. Marques, Physica A 172, 
285 (1991). 

E. E. Dormidontova, Macromolecules 32, 7630 (1999). 
C. N. Likos, H. Lowen, M. Watzlawek, B. Abbas, O. 
Jucknischke, J. Allgaier and D. Richter, Phys. Rev. Lett. 
80, 4450 (1998). 

E. g. M. Born and K. Huang, Dynamical Theory of Crys- 
tal Lattices (Oxford University Press: Oxford, 1954). 

F. A. Lindemann, Z. Phys. 11, 609 (1910). 

J. P. Hansen and L. Verlet, Phys. Rev. 184, 151 (1969). 
T. A. Witten, P. A. Pincus and M. E. Cates, Euro. Phys. 
Lett. 2, 137 (1986). 

B. J. Alder, W. G. Hoover and D. A. Young, J. Chem. 
Phys. 49, 3688 (1968). 

P. Ziherl and R. D. Kamien, Phys. Rev. Lett, 85, 3528 
(2000). 

G. M. Grason, Phys. Rep. 433, 1 (2006). 

L. Leibler, H. Orland and J. C. Wheeler, J. Chem. Phys. 
79, 3550 (1983). 

M. Imai, I. Yoshida, T. Iwaki and K. Nakaya, J. Chem. 
Phys. 122, 044906 (2005). 

E.g. J. N. Israelachvili, Intermolecular and Surface 
Forces, 2nd ed. (Academic Press: London, 1992). 



